Tunable Holstein model with cold polar molecules 
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We show that an ensemble of polar molecules trapped in an optical lattice can be considered as 
a controllable open quantum system. The coupling between collective rotational excitations and 
the motion of the molecules in the lattice potential can be controlled by varying the strength and 
orientation of an external DC electric field as well as the intensity of the trapping laser. The system 
can be described by a generalized Holstein Hamiltonian with tunable parameters and can be used 
as a quantum simulator of excitation energy transfer and polaron phenomena. We show that the 
character of excitation energy transfer can be modified by tuning experimental parameters. 



Introduction.- Many important biological processes 
involve energy transfer between complex molecules in 
mesoscopic aggregates [H, Energy can generally be 
transferred incoherently via direct two-molecule interac- 
tions or through the emergence of collective coherence 
described by excitons [3[. Incoherent energy transfer re- 
sults from the interaction of excitons with phonons. It is 
a highly debated open question whether exciton-phonon 
interactions conspire to ensure the most efficient and uni- 
directional energy transfer in biological systems 
The exciton-phonon interactions can be described by a 
Holstein model @, [1] ■ It is difficult to calculate numer- 
ically the full energy spectrum for this polaron model 
in a complete range of interaction parameters Q . There- 
fore, it is necessary to design an experimentally accessible 
many-body quantum system that would be described by 
the Holstein Hamiltonian with tunable parameters and 
arbitrary dimensionality. Such a system could be used for 
quantum simulation of excitation energy transfer (EET) 
in complex molecular aggregates and polaron phenomena 
in general. 

There is currently growing interest in using ultracold 
atoms trapped on an optical lattice for quantum sim- 
ulation of condensed matter physics [13, [HI- Ultra- 
cold atoms offer the possibility of designing systems that 
are well described by model Hamiltonians such as the 
Bose-Hubbard Hamiltonian 12]. It was recently shown 
that polar molecules trapped on an optical lattice pro- 
vide new possibilities for quantum simulation due to the 
presence of long-range dipole-dipole interactions [lljl - iis] . 
Here, we consider the interaction of rotational excitons 
with phonons in an ensemble of ultracold polar molecules 
trapped in an optical lattice. We show that, although the 
translational motion of molecules is largely determined 
by the intensity of the trapping laser, the dipole-dipole 
interaction can be made large enough to couple the dy- 
namics of rotational excitons with the lattice vibrations. 
We demonstrate that by tuning the trapping laser in- 
tensity and an applied DC electric field, the strength of 
the exciton-phonon coupling can be controlled, and that 
the character of EET can be modified dynamically from 



coherent to incoherent. 

Molecular crystal Hamiltonian - We consider an array 
of polar molecules in the rovibrational ground state, 
trapped in a 3D optical lattice with one molecule per 
lattice site and no tunneling of molecules between sites 
pH [l2| . Trapping of 1 £ diatomic molecules on an optical 
lattice has recently been demonstrated [l6l - [l8| . For the 
lowest bands of the periodic lattice potential, molecules 
vibrate harmonically around the equilibrium positions Rj 

[Eta. 

The trapping strength of the optical lattice is one ex- 
perimental parameter to control the system. Another 
parameter can be introduced by applying a DC electric 
field E. In a weak DC field, the rotational ground state 
\g) w a\N = 0,Mjv = 0) +b\N = l,M N = 0) and the ex- 
cited state |e) « b\N = 0,M N = 0) - a\N = 1,M N = 0) 
constitute an isolated two-level system [3, [H[ . The field- 
free rotational states \N, M^) are eigenstates of the rigid 
rotor Hamiltonian Hr = B e N 2 , where B e is the rota- 
tional constant. The states \g) and |e) are eigenstates of 
the Hamiltonian H-qq — H-r ~ d-E, where d is the electric 
dipole operator. The coefficients a and b are functions of 
the DC field strength E. The electric dipole-dipole inter- 
action Vi(ri, Yj) couples the rotational states of molecules 
in different lattice sites. For lattice site separations ol ~ 
500 nm and molecules with a permanent dipole moment 
d > 1 Debye, the characteristic energy of the dipole- 
dipole interaction Vdd = d 2 j a\ is a few tens of kHz. In 
the two-molecule subspace S = {\g,g}, \g, e), |e, g)}, the 
dipole-dipole operator Vi(ri,Tj) has the following matrix 
(9i,9j\Vi\9i,gj), V*j 9 = (ei,gj\Vi\ei,gj), 
These integrals can be eval- 



can be written as V, 
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and Jij = (gi,ej 
uated as in Ref. 

£V|r*-r,| 3 . 

An ensemble of polar molecules in an optical lattice 
can be represented by a Hamiltonian of the form % = 
H p h + H cx + Hi n t ■ The first term, 
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describes phonons associated with the oscillatory mo- 
tion of the molecules in the lattice potential V g (ri,rj) — 
moj 2 (ri — Ri) 2 /2 + U g /\ri — rj\ 3 . The first contribution 
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to V g depends on the intensity of the trap ping laser that 
determines the trapping frequency ujq [10| and the molec- 
ular mass m. The second term in Eq. (JT]) depends on 
the strength of the dipole-dipole interaction and couples 
the motion of molecules in different sites. The opera- 
tor a k v creates a phonon in mode k with polarization 
v = x,y,z. A competition betwen the laser trapping 
force /i oc mwj and the dipole-dipole force fdd oc U g /a 5 L 
determines the phonon frequency Uk.w For example, for 
a homogeneous ID array uik = woyl + ^PtPOj where 
P = (Ug/a 5 L )/mu$ and7(fc) = Ei>o[ 1_cos 0' fc )]/i 5 ' with 
k within the first Brillouin zone. For finite arrays, the 
phonon frequencies are obtained from the eigenvalues of 
the force constant matrix T given by the Hessian of the 
lattice potential V = J2i j ^s( r i7 r j ) 20]. The phonon 
spectrum in optical lattices is gapped for any value of p, 
i.e., LJk — > wo as k — > 0, which resembles optical phonons 
in solids. In the limit p <C 1, the spectrum is disper- 
sionless (Einstein oscillators) . This limit can be achieved 
either by increasing the lattice depth or by decreasing 
the DC field strength. In experiments, the Gaussian pro- 
file of the trapping beams usually generates an additional 
global harmonic potential Vh(rj) = mo;^r|/2. This po- 
tential can be included in V g (Ti,Tj). However, for the 
optical lattice potentials considered here Vh can be ne- 
glected since Wh/^o < 10~ 2 10]. 

The second term in the total Hamiltonian TL describes 
collective rotational excitations (excitons) and is given 

by his p 

4x = $>e 9 + Di)B\Bi + £ J id B\Bj ■ (2) 

i 

In the present work, we neglect non-linear exciton- 
exciton interactions, which is a good approximation for 
a small number of rotational excitations. The transi- 
tion operator B\ = \ei)(gi\ creates a rotational excita- 
tion \g) —> e) in site i. The first term in Eq. ([2]) con- 
tains the excitation energy at each site, which is a sum of 
the single-molecule rotational splitting e eg and the site- 
dependent shift Di = J2j^i Dij due to the dipole-dipole 
interaction between molecules, where = V^ 9 — V? 3 . 
The second term in Eq. (£Z$ describes the hopping of the 
rotational excitation between sites. The integrals Dij 
and Jij in Eq. ([2} are evaluated for molecules fixed at 
their equilibrium positions R^. 

The final contribution to the system Hamiltonian % is 
the exciton-phonon interaction. This effective coupling 
between internal and translational degrees of freedom of 
the molecules can be written as 

Hy A =Y,9D i {al v + a kv )B\B i 

k,i 

k Lj^i 

which derives from the Taylor expansion of the integrals 
Dij and Jij in Eq. ^ up to linear order with respect 
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FIG. 1. (color online) Nearest-neighbour couplings J12 (solid 
line), D12 (dashed line), and (dot-dashed line): (a) de- 
pendence on the DC field strength parameter A = dE/B e for 
a DC field E = Ez perpendicular to the one-dimensional ar- 
ray; (b) dependence on the angle 9 between the electric field 
and the array, for A = 1. Energy in units of Vdd = d 2 /a\, 
where d is the permanent dipole moment and ol is the lattice 
constant. B e is the rotational constant of the molecule. 

to the small variation of the relative distance 8t^ from 
its equilibrium value, due to the motion of the molecules 
in their local potentials. The exciton-phonon coupling 
parameters gff. and gj v . can be obtained analytically 
from the gradient of the dipole-dipole potential and de- 
pend on the phonon mode (k, v). For finite arrays, the 
mode dependence can be obtained from the eigenvec- 
tors of the force constant matrix T . The first term in 
Eq. © represents phonon-modulated site energies, pro- 
portional to Dijd^ 1 , and the second term corresponds 
to phonon-modulated hopping of an excitation, propor- 
tional to JijCtjj ■ Although the energy shift is a small 
perturbation to the rotational spectrum, i.e., <C e e g, 
its fluctuation with the motion of the molecules in their 
local potentials can lead to dynamical localization of ex- 
citons [1,0. 

The system Hamiltonian H in Eqs. CD)-© constitutes 
a generalized polaron model. In the limit g^ ^> g k f._ , 
it corresponds to the standard Holstein model @, 
extensively used to study energy transfer in molecular 
crystals and photosynthetic complexes 0, @, polaron 
physics in solids (22[, and, since recently, quantum infor- 
mation processing in dipolar gases [2lJ. In the opposite 
limit g^[ <C g k ] V . , the Hamiltonian corresponds to the 
Su-Schricffcr-Hceger (SSH) model of particle-boson cou- 
pling, introduced to describe electrons in one-dimensional 
chains of polyacetylene [23|, H3] ■ It is easy to include ad- 
ditional effects such as the anharmonicity of the optical 
lattice potential [25[, quadratic exciton-phonon coupling 
[2(| and exciton-exciton interactions Q in the Hamilto- 
nian. 

A similar polaron model T-L was used in Ref. [2l[ to 
describe 2D self-assembled crystals of polar molecules in 
the context of quantum information processing. Quasi- 
2D crystals are predicted to form at temperatures T ~ 10 
nK, when the molecular kinetic energy is smaller than 
Vdd- A strong transverse confinement and a DC electric 
field perpendicular to the crystal plane are needed to sta- 
bilize the crystal against attractive dipole-dipole interac- 
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FIG. 2. (color online) Polaron shift A_B 3 as a function of the 
trapping frequency u)o for an array of 10 LiCs molecules sep- 
arated by 400 nm. Curves are shown for electric fields of 9 
kV/cm (solid line), 2 kV/cm (dashed line), and 0.6 kV/cm 
(dotted line). Panels (a) and (b) correspond to a field per- 
pendicular and parallel to the array, respectively. In panel 
(a) we have \AE g \ > J 12 = 6.73 kHz, for E = 9 kV/cm 
and uso/2n < 20 kHz, which is a signature of strong exciton- 
phonon coupling. 

tions. Self-assembled crystals exhibit acoustic phonons 
and the exciton-phonon interaction is always strong due 
to the presence of low frequency modes. These con- 
straints limit the range of the Hamiltonian parameters 
(J, D, g'l), gj) that can be explored with such systems. 
In the system proposed here, molecules are stabilized 
against collisional losses by the optical lattice potential 
for any orientation of the DC electric field, and for any 
dimensionality. In addition, the coupling parameters can 
be tuned in a much wider range of values, as demon- 
strated below. 

Tunable exciton-phonon coupling.- We now specialize 
our discussion to a finite ID array of polar molecules. 
Effective lower-dimensional arrays in a 3D optical lattice 
can be generated when the dipole-dipole interaction is 
significant along one or two axes of the lattice only. We 
consider the interaction of a single exciton with harmonic 
phonons. 

The DC electric field modifies the single-molecule 
states \g) and |e), and therefore the value of the dipole- 
dipole couplings J12, , and Vf 2 9 [HI- This is shown in 
Fig. [Ha) for a DC field perpendicular to the axis of the 
array. The magnitude of J 12 decreases with increasing 
field strength, whereas D\i and Vf 2 9 increase. For large 
DC fields, Eq. (|3]) reduces to the Holstein polaron model 
plus a small correction due to the finite value of gj 12 . In 
the limit of weak DC fields, D12 and Vf 2 9 are vanishingly 
small due to parity selection rules, and Eq. ([3]) reduces to 
the SSH polaron model with Einstein phonons. The val- 
ues of the dipole-dipole matrix elements also depend on 
the angle 9 between the DC electric field and the molecu- 
lar array. For our chosen rotational subspace with projec- 
tion Mjv = along the electric field axis, the parameters 
-D12, J12, and Vf 2 9 are proportional to (3 cos 2 9 — 1) and 
vanish for 9 ps 54.7°, as shown in Fig. [TJb). 

The strength of the exciton-phonon coupling can be 
controlled by tuning the trapping laser intensity and the 
DC electric field. The coupling constants in Eq. © can 



be written as guj oc yj (1 /muJo) (Ax2 / a) , for A = D, J. In 
order to quantify the strength of this coupling for realistic 
trapping conditions, we analyze the eigenvalues of the 
Holstein Hamiltonian for a given molecular species in a 
finite ID array. We diagonalize the total Hamiltonian 
% numerically for an array of a molecules in the site 
basis \gx, ■ ■ ■ , e», . . . ,g a )\vi, v%, . . . , v a ), where Vk is the 
occupation number of the phonon mode k. The phonon 
basis is truncated by including states with up to a given 
phonon occupation ^ max - The value of i/ max is increased 
iteratively until the calculated observable is converged. 
We partition the Hamiltonian as T~L = Hq + Hi nt , where 
H n = H ex + H p h. The ground state energy E g of the non- 
interacting Hamiltonian Hq is chosen as a reference. Any 
interaction between rotational excitons and the lattice 
vibrations of the molecules shifts the ground state E g 
towards lower energies. As an illustrative example, we 
show in Fig. [2] the shift AEJo for a finite array of LiCs 
molecules (d = 5.5 Debye [27|) separated by ai = 400 
nm, as a function of the optical lattice trap frequency 
loq. Each curve corresponds to a different field strength 
E < 10 kV/cm. For small trapping frequencies (weaker 
lattices), we find AE g > J12, which is a signature of 
strong coupling [8( . The strength of the exciton-phonon 
coupling is larger for a DC field parallel to the array than 
for any other field orientation. 

Application to energy transfer in a phonon bath. - One 
of the possible applications of a tunable Holstein Hamil- 
tonian using cold polar molecules is the simulation of 
excitation energy transfer processes (EET) that occur 
in molecular crystals and light-harvesting complexes at 
room temperature P, H, [H, H H^. For example, let us 
consider two LiCs molecules in an optical lattice with 
ol — 400 nm and trapping frequency uj /2it = 10 kHz, in 
a DC field of 10 kV/cm perpendicular to the intermolec- 
ular axis. The two normal modes of lattice vibration 
have frequencies io\ — loq and C02 = 2.8ujq. The lower fre- 
quency mode does not couple to excitons because it does 
not change the relative distance between molecules. The 
Hamiltonian parameters for the higher frequency mode 
are gD 12 /h — 12.7 kHz and gj 12 /h = —2.33 kHz. This 
gives the ratio <7d 12 /Ji2 ~ 2, which can also be found 
in the Fenna-Mathews-Olson photosynthetic complex, 
where electronic excitations are believed to be locally 
coupled to phonons at each site For a weak DC field 

E = 0.5 kV/cm, the first term in Eq. © is negligible, 
and non-diagonal coupling dominates {gj 12 / J12 ~ 0.6, for 
u>o/2tt = 10 kHz), which may allow for tests of the role 
of spatial non-local phonon correlations in the dynamics 
of EET 

In order to model rotational EET using the generalized 
polaron Hamiltonian in Eqs. ([TJ-©! we define an ini- 
tial wavefunction |\&(0)) describing the coupled exciton- 
phonon system, using the product basis described above. 
We propagate the time-dependent Schrodinger equation, 
using the Hamiltonian T-L, and construct the total den- 
sity matrix p(t) — (t)} (t)\ at each time step. We 
then obtain the reduced density matrix in the exciton 
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FIG. 3. (color online) Excitation energy transfer in an array 
of five LiCs molecules in a DC electric field perpendicular to 
the array: (a) Evolution of the excitation probability pi(t), 
when no phonons are present; (b) The same as in (a), but 
with phonons in an optical lattice with trapping frequency 
va = oJo/2tt varying in time as indicated in the inset. The 
field strength is 0.5 kV/cm. 

subspace pE(t) = Tr v ;b{/5(i)} by tracing over the states 
in the truncated phonon basis. The diagonal elements 
of the reduced density matrix are the time-dependent 
probabilities Pi(t) for a molecule in site i to be in the 
rotational excited state | e^) . We apply this procedure to 
a ID array of five LiCs molecules separated by 400 nm, 
in a weak DC field E — 0.5 kV/cm perpendicular to the 
array. The rotational excitation is initially in molecule 
1. Single-site excitation and probing of rotational states 
can be achieved by app lying an electric field gradient as 
described in Rcf. [30]. The phonon bath has initially 
zero temperature. In Fig. [31 we show the dynamics of 
the excitation at site 1 for different trapping frequen- 
cies u!q/2tt, keeping E constant. In Fig. EJa) we set 
loq — > oo, thereby neglecting the coupling to phonons. In 
this limit, the excitation energy oscillates coherently be- 
tween the molecules of the array [l| , with a transfer time 
between adjacent sites t\i — h/\Ji 2 \ ~ 10 [is. When 
the exciton-phonon coupling is turned on, by decreasing 
the trapping frequency in this case, the transport of the 
excitation is suppressed due to a competition between 
the exciton-phonon interactions and the excitonic energy 
transfer. When the exciton-phonon coupling is strong 
enough, all the site populations Pi(t) approach an equi- 
librium value. This is shown in Fig. [3jb), where only 
the population in site 1 is presented for simplicity. We 
obtain a similar behaviour when ujq is kept constant and 
the DC field strength is dynamically tuned from weak to 
strong. 



Discussion.- We have shown that the translational and 
rotational states of polar molecules trapped in an optical 
lattice can interact in the presence of a DC electric field. 
This interaction is described by a generalized polaron 
model with tunable parameters. All the terms in the 
Hamiltonian can be dynamically tuned by varying the 
DC electric field or the lattice laser intensity. 

These results suggest the possibility of using cold po- 
lar molecules for quantum simulation of excitation en- 
ergy transfer processes that occur at high temperatures in 
solids and mesoscopic systems of biological interest PHI] . 
Despite the physico-chemical differences between these 
systems and cold polar molecules, they are described by 
the same Hamiltonian, and it should be possible to use 
the latter to experimentally study the role of environmen- 
tal noise in the efficiency of quantum transport [5l l3l[l32l ]. 
In addition, polar molecules can also be used to explore 
the strong coupling regime of exciton-phonon interaction, 
which may provide insight into some of the open ques- 
tions of polaron physics, such as the role of the lattice 
dynamics in the mechanism of high-T c superconductivity 

The interactions of atomic spin states with phonons 
are used in experiments with ultracold trapped ions to 
generate phonon-mediated gates for quantum comput- 
ing [33j]. Similarly, the interactions of rotational states 
with lattice vibrations described in this work can be used 
to produce novel quantum gates. The detrimental ef- 
fects of decoherence can be minimized by increasing the 
trapping frequency. Phonon-mediated interactions may 
also introduce new adjustable parameters to the spin- 
lattice Hamiltonians that can be simulated with ultra- 
cold molecules [H, |34j. In particular, the controllable 
coupling to phonons can be exploited to study the effects 
of dynamically tunable decoherence on topologically or- 
dered states that can be created with ultracold molecules 
13]. In summary, the proposed system may provide an 
experimental tool to explore the dynamics of quantum 
networks whose essential features are replicated in sev- 
eral areas of condensed matter physics. The ability to 
tune the system-environment coupling dynamically may 
also stimulate new studies of the dynamics of open quan- 
tum systems. 

Note.- After submission of this work, quantum simula- 
tion of EET using superconducting qubits was proposed 
in Ref. [35|. 
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